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ABSTRACT 

We investigate the importance of interactions between dark matter substructures 
for the mass loss they suffer whilst orbiting within a sample of high resolution galaxy 
cluster mass Cold Dark Matter haloes formed in cosmological N-body simulations. We 
have defined a quantitative measure that gauges the degree to which interactions are 
responsible for mass loss from substructures. This measure indicates that interactions 
are more prominent in younger systems when compared to older more relaxed systems. 
We show that this is due to the increased number of encounters a satellite experiences 
and a higher mass fraction in satellites. This is in spite of the uniformity in the 
distributions of relative distances and velocities of encounters between substructures 
within the different host systems in our sample. 

Using a simple model to relate the net force felt by a single satellite to the mass 
loss it suffers, we show that interactions with other satellites account for ~ 30% of 
the total mass loss experienced over its lifetime. The relation between the age of the 
host and the importance of interactions increases the scatter about this mean value 
from ~ 25% for the oldest to ~ 45% for the youngest system we have studied. We 
conclude that satellite interactions play a vital role in the evolution of substructure in 
dark matter halos and that a significant fraction of the tidally stripped material can 
be attributed to these interactions. 

Key words: galaxies: clusters - galaxies: formation - galaxies: evolution - n-body 
simulations 



1 INTRODUCTION 

It has been understood for some time that the structure 
of a galaxy can be affe cted by tidal interaction s with its 
close neig hbour(s) (e.g. iToomre fc Toomrelll972ft : tell-tale 
signs such as tidal tails and disturbed morphologies pro- 
vide a visible record of these encounters. Around our own 
Galaxy, there is substantial evidence for its tidal interac- 
tion with the Small and Large Magellanic Clouds (SMC and 
LMC ), the consequences of which have been studied in detail 
(e.g. iLin et all Il995t lOh et alJ Il995t iGardiner fc Noguchil 
19961: lYoshizawa fc Noguchil 120031: iBekki fc Chibal 120051: 
Mastropietro et alJ 120051 : IConnors et alJ l2005h . Further- 
more, an increasing number of studies have uncovered evi- 
dence for tidal stripping - in the form of stellar streams - 
in the Galactic halo (e.g. iHelmil 12004 llbata et al-lEioll) : 
these streams represent material that has been stripped 
from infalling satellites as they are disrupted by our Galaxy. 



The detection of such streams will become more common- 
place in the coming years as the sen s itivity of surveys im - 
prove (e.g. lOdenkirchen et al-lEoo! iNavarro et al]l2004l1 . 
but there are already exam ples of stellar stre ams further 
afield, such as around M31 (llbata et alj|2004) . Moreover, 
iMihos et alJ (120051) recently reported the discovery of in- 
tracluster light in the Virgo cluster, revealing several long 
(>100 kpc) tidal streamers. 

These results represent compelling evidence that 
satellite galaxies tidally interact with their more massive 
hosts, and consequently lose some fraction of their mass. 
The effect of a satellite's interaction with its host and the 
mass loss it suffers has been studied in some detail (e.g. 
lHavashi et al J 120031. and it can be argued that it is rela- 
tively well understood. In comparison, the importance of a 
satellite galaxy's interactions with other satellite galaxies, 
the nature of these interactions and the contribution they 
make to its mass loss is less well understood. There is 
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evidence to sugg est th a t tidal intera ctions between satellite 
galaxies occur; IZhaol il998f) and llbata fe Lewis! (Il998h 
investigated whether the Sagittarius Dwarf galaxy could 
have experienced an encounter with the SMC and LMC 
some 2-3 Gyrs ago, while the disturbed HI distribution 
noted by lYun et all l)l994h in the M81 group is highly 
sugge s tive o f tidal interactions between the group galaxies. 
iGotol i2005T) argues that tidal interaction between galaxies 
is the dominant mechanism driving cluster galaxy evo- 
lution and underpins the Butcher-Oemler effect and the 
morphology-density relation. 

It has been understood for some time that dark matter 
haloes must play an important dynamical role in encounters 
between galaxi es because th ey significantly reduce the merg- 
ing timescale l)Barneslll98Sf l. Examples of tidal interaction 
and merging are observed in relatively low-density environ- 
ments (i.e. the field), but how reasonable is it to expect 
that interactions should be more common in higher density 
environments such as galaxy groups and clusters? Tidal in- 
teractions have been proposed as a mechanism for galaxy 
transformation in gal axy clusters, such a s the "harassment" 
scenario envisaged bv lMoore et alJ <ll998l) , but what does the 
favoured paradigm for cosmological structure formation, the 
Cold Dark Matter model, predict? 

The aim of this paper is to quantify the importance of 
satellite-satellite encounters and to assess their impact upon 
the transformation and mass loss of the substructure popu- 
lation within the context of the Cold Dark Matter (CDM) 
model. We have drawn on a sample of high resolution cos- 
mological N-Body simulations of cluster mass dark matter 
haloes and analysed the interactions of the substructure 
haloes (hereafter subhaloes) both with the host halo and 
with other subhaloes. We associate these subhaloes with the 
hosts of satellite galaxies (but see iGao et al ] |2004T) and in 
what follows we use the terms subhalo and satellite (galaxy) 
interchangeably. The fine time sampling of our simulations 
allow us to follow the time evolution of subhalo properties in 
detail, and so we can determine the relative contributions of 
the host and the other subhaloes to changes in a subhalo's 
structure. 

In a previous study llKnebe et al] l2004) we quanti- 
fied the frequency of encounters between subhaloes orbit- 
ing within a common CDM cluster mass halo, considering 
the period between the halo's formation redshift 1 and the 
present day. We found that, on average, 30% of the "satellite 
galaxy" population experienced at least one encounter per 
orbit with another satellite galaxy. This result was sensitive 
to the age of the host halo, with a clear trend for more en- 
counters in younger systems. We also reported a correlation 
between the number of encounters and halocentric radius - 
satellite galaxies closer to the centre of the host were mea- 
sured to experience more interactions, although we note that 
this simply reflects the increasing spatial density of satellites 
with decreasing radius within a host halo. 

The principal shortcoming of the approach adopted in 
iKnebe et al] i2004h is that we neglected the relative veloc- 



We defined this to be the redshift at which the mass of the 
most massive progenitor was half the system's present day mass; 
this was typically z ~ 0.5 for the haloes we examined. 



ities of the satellite galaxies; our satellites may have expe- 
rienced encounters, but we had no information about their 
specific nature, i.e. were they fast or slow? Such information 
is important when considering the impact on the satellite's 
structure. In this present study, we elaborate on that work 
by including information about the relative velocities of the 
satellites. In other words, we can now estimate the impor- 
tance of encounters in addition to the frequency with which 
they occur, allowing us to differentiate between slow en- 
counters, which we expect to be extremely disruptive to the 
satellite structure, and fast encounters, whose impact are 
likely to be minimal. We define a quantitative measure for 
interactions, which we call the integral interaction measure, 
based upon the force acting on a satellite over a given period 
of time, i.e. the (induced) momentum change. Whereas be- 
fore we could examine the number of encounters a satellite 
galaxy experienced per orbit, we may now study how the 
instantaneous force due to encounters acting on a satellite 
galaxy varies along its orbit and how this correlates with 
mass loss, thus providing a natural measure of the impor- 
tance of mutual interactions between satellite galaxies. 

In what follows, we motivate our choice of the inte- 
gral interaction measure (hereafter IIM) as a gauge for 
the importance of interactions between satellite galaxies, 
and present the results of our analysis of a suite of high 
resolution cluster mass haloes that formed assuming the 
ACDM cosmology. We demonstrate the suitability of the 
IIM for our purposes by performing a series of experiments 
with "cleaned" simulations, in which we track the detailed 
mass loss history of a single satellite galaxy in a host halo in 
which the substructure has been removed. Finally, we com- 
pare and contrast our results with those of previous studies, 
and comment on their observable consequences. 



2 THE SIMULATIONS 

Our analysis is base d on a suite of eight high-resolution 
A-body simulations llGill et ail l2004al lbfl carried out using 
the publicly available adaptive mesh refinement code MLAPM 
jKnebe. Green fc BinnevlboOll) in a standard ACDM cos- 
mology (fi = 0.3, n x = 0.7, Q b h 2 = 0.04, h = 0.7, <r 8 = 0.9). 
Each run focuses on the formation and evolution of a dark 
matter galaxy cluster containing of order one million parti- 
cles, with mass resolution 1.6 x 10 8 h^ 1 Mq and force reso- 
lution ~2/i~ 1 kpc which is of the order 0.05% of the host's 
virial radius. These simulations have the required resolution 
to follow the satellites within the very central regions of the 
host potential (5=5—10% of the virial radius) and the time 
resolution to resolve the satellite dynamics with good ac- 
curacy (At ~170 Myrs). Such temporal resolution provides 
of order 10-20 timesteps per orbit per satellite galaxy, thus 
allowing these sim ulations to be used in a previous paper 
jGill et aljl2004bh to accurately measure the orbital param- 
eters of each individual satellite galaxy. 

Substructure within these haloes is identified using the 
halo finder MHF (MLAPM's-halo-finder) . MHF is based upon the 
A^-body code MLAPM and acts with exactly the same accuracy 
as the A-body code itself; it is therefore free of any bias 
and spurious mismatch between simulation data and halo 
finding precision arising from numerical effects. We applied 
MHF to each of our eight host halos at their formation time 
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Figure 1. The integral interaction measure as a function of satel- 
lite mass. 



which is the redshift zf orm where the halo contains half of 
its present day mass. We track the orbits of each of the 
satellites identified within and around the host halo from 
2form until z = and follow the evolution of their properties 
in great detail. For further details relating to the properties 
of the s ate llite galaxies, we refer the rea der to the lGill et alJ 
i2004eJbh : iGill. Knebe fc Gibson! J2005h series of papers. 



3 THE ANALYSIS 

In what follows, we have considered only those satellites 
that have completed at least one full orbit within their host 
halo, corresponding to more than 70% of the subhaloes. 
The number distribution of orbits peaks at about 1-2 or- 
bits with the tail extending to as many as four orbits for the 
older host halos. A detailed discussion of the orbital prop- 
erties (amo ngst others) of the substructure population can 
be found in lGill et alJ J2004al) . 

We restricted our sample of satellites to those that con- 
tain at least 100 particles, which translates to a minimum 
mass of M sat ^ 2x 10 10 h' 1 Mq. To ensure that our results 
are not affected by resolution effects, we checked that all 
results presented below are recovered when the lower mass 
limit is gradually increased; that is, we considered additional 
lower mass cuts corresponding to 200 and 400 particles and 
we can confirm that our results are unaffected. 



3.1 Integral Interaction Measure 

We begin by calculating the forces acting on each satellite 
galaxy i for each available snapshot of the simulation, treat- 
ing it as a point particle with mass m. The force F^ ost ex- 
erted by the host halo and the force Fg at exerted by all other 
satellites are given as follow 
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where rrn is the mass of satellite i and Mh os t(< n) the mass 
of the host interior to the satellite distance Ti. We need to 
stress that both these formulae assume spherical symmetry 
and hence are only approximations to the "true" forces. 
We define a so-called (dimensionless) "integral-interaction- 
measure" - IIM - for each individual satellite galaxy as fol- 
lows: 
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where we integrate over a time interval [0,T], which is the 
time satellite i has spent within its host's virial radius. Here 
we also note that due to our definition the IIM values scale 
linearly with the "average satellite mass". The discrete na- 
ture of the time sampling of our data requires that the in- 
tegral should be expressed as the following summation: 



IIM 1 



tnow 



F sat (t m ) 

Fhost(tm) 



At 



(3) 



where t n0 w is the age of the Universe at redshift z = 0, ti the 
age of the Universe when the satellite enters the host halo, 
and At the time difference between two consecutive outputs. 
We average the forces exerted by both the other satellites 
and the host halo over the consecutive outputs, i.e. [t — 
At/2, t + At/2], or "mid-point integration" of equation @: 
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The integral-interaction-measure equation © can now be 
used as a quantitative measure for the relative strength of 
satellite-satellite encounters. 



3.1.1 Application of the Integral Interaction Measure 

In Fig. Q we present the integral-interaction-measure IIM, 
as defined by equation 10, for each satellite in our suite of 
eight host halos plotted as a function of satellite mass. This 
figure suggests that there is no clear trend for interactions 
to correlate with mass, as we might have expected; it would 
be rather surprising to find that, for instance, high-mass 
satellites tend to interact more prominently than low-mass 
ones (or vice versa). 

The most striking feature of Fig. is the apparent rise 
of the IIM values as a function of decreasing age for the host 
halos: the halos are ordered in age with halo #1 being 8.3 
Gyrs old and halo #8 a mere 3.4 Gyrs. 
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Table 1. Best fit parameters of IIM distribution to lognormal. 



halo 


age [Gyr] 


HM peak 


IIM 


fo 


# 1 


8.30 


0.067 


0.118 


0.756 


#2 


7.55 


0.137 


0.199 


0.614 


#3 


7.16 


0.126 


0.227 


0.767 


#4 


7.07 


0.162 


0.270 


0.716 


#5 


6.01 


0.287 


0.464 


0.692 


#6 


6.01 


0.221 


0.307 


0.575 


#7 


4.52 


0.535 


0.789 


0.623 


#8 


3.42 


0.672 


1.021 


0.646 



This can be better viewed in Fig. [5] where we plot the 
distributions of the integral interaction measure IIM. For all 
our eight host halos these distributions have been fitted with 
a log-normal function 



n(IIM) 



V 27r °o 



exp 



ln 2 (IIM/IIM ) 



2al 



(5) 



corresponding best fit parameters along with the halo age 
are listed in TableQwhere IIM pca k = IIMo exp(— a 2 ) for a 
log-normal distribution. 

The increase of the IIM with decreasing a ge of the host 
is con sistent with the behaviour observed in iKnebe et alJ 
(2004), in which it was noted that the tail of the distribu- 
tion of the number of encounters per orbit extended to larger 



values for younger host systems. However, the r esult implied 
by the IIM is distinct from that presented in lKnebe et alJ 
(2004) in the sense that we are considering the net force 
acting on a subhalo over some time interval, whereas we 
previously considered encounters as events in which a pair 
of subhaloes were spatially coincident. This raises the ques- 
tion of whether or not the IIM is a reasonable measure of 
interactions, and in particular, if it could simply be the case 
that it is dominated by single encounters. 

We investigate this in Fig. |3 where we examine the 
correlation between the IIM and the number N cnc of "tidal 
encounters" as quantified by calculating the tidal radius 
of a given satellite induced by one of the other satellites 
jKnebe et al.l l2004) . Whenever the tidal radius becomes 
smaller than the radius 2 of the satellite we increment a 
counter iV cnc for that particular satellite that keeps track 
of the number of (perturbing) interactions with companion 
satellite galaxies. 

Fig. El clearly indicates that there is little (if any) cor- 
relation between the number of satellite-satellite encounters 



2 We define the radius of a satellite either to be the virial ra- 
dius, i.e. the radius where the mean averaged density (measured 
in terms of the cosmological background density p;,) drops below 
A v i r (,z), or the truncation radius, i.e. the point where the satel- 
lite's density profile rises again due to the embedding in the host 
halo. 
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and the integral interaction measure for a single satellite. 
Moreover, also the scatter about the mean IIM value in each 
iVcnc bin is not affected by the actual number of encounters 
experienced by the satellite. This strongly suggests that the 
IIM value is not dominated by single events but rather is 
a cumulative quantity that is accrued over the lifetime of a 
satellite. However, we stress that there is a correlation be- 
tween the w idth of the distribu tion of encounters per orbit 
(c.f. Fig.3 in lKnebe et"aill200^ and the peak IIM value; al- 
though the IIM is not driven by single violent encounters, 
the greater the number of such events, the higher the IIM 
of the satellite. However, these arguments are based upon 
the assumption that the "strength" of individual encounters 
is more or less equal. It still appears possible for one single 
strong encounter to dominate the value of IIM. 

Another factor possibly affecting the observed correla- 
tion between host halo age and interaction measure IIM is 
the mass fraction of satellites. A simple check indicates that 
the younger the host the higher the fraction of mass locked- 
up in satellite galaxies. This suggests that the IIM values 
may in fact be influenced by the most massive subhaloes. 
We will come back to this point later in Section l3.4l but can 
already confirm that the distributions presented in Fig. [5] 
practically remain unaltered if we discard all satellites less 
massive than 1% of the host's virial mass 3 , denoting the 
importance of massive subsystems. 

In addition, we have investigated whether or not there 
exists a relation between then IIM and either the eccentric- 
ity of a satellite's orbit or its pericentric distance, but we 
do not find strong evidence for such a correlation. Although 
we observed a significant drop in the number of encounters 
per orbit with increasing distance from the host's centre, we 
find no comparable result for IIM. This indicates that satel- 
lites "encounter" each other with greater frequency closer to 
the centre of the host, but that such encounters occur with 
high relative velocities and so cause little structural dam- 
age. Encounters in the central regions are therefore no more 
damaging than those in the outer regions. We elaborate upon 
this in greater detail in the following section. 

3.2 Distributions of Relative Encounters 

We have calculated the distribution of satellite-satellite dis- 
tances Di,j as well as the relative speed of satellite pairs Vi l j l 
for all available outputs in-between formation redshift 2f orm 
of the host and z = 0, and show the resulting distributions 
of in Fig. 2] and Fig. |5] respectively. 

In Fig. |3] we have plotted the relative separation Dij 
normalised by the sum of the two virial radii of the respec- 
tive satellites, i.e. Ri and Rj; a value of (Ri + Rj)/Dij > 1 
corresponds to a distance of the satellites-satellite., pair for 
which the "virial spheres" of the satellites i and j are over- 
lapping. We note that the distributions can be fitted by a 
log-normal distribution 

3 One needs to bear in mind that the mass spectrum of sub- 
haloes extendsdown to as low as 10 — 4 — 10 — 5 M v i r host (e.g. 
iDeLucia et al.ll2004) 
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Figure 4. Distribution of the distance of two satellites normalized 
by the sum of their individual radii; data from each host halo are 
stacked for all available outputs. 



where x = (Ri + Rj) / Dij . Fig.^Jis accompanied by Tabled 
where we summarize the best-fit parameters. Despite the 
age-IIM relation found in the previous Section 13.11 we do 
not observe any trend fo r relative distances t o increase (or 
decrease) with halo age. iTormen et al.l il998t) performed a 
similar analysis, but their respective distance distribution 
peaks for values corresponding to distances smaller than the 
sum of the two individual radii indicating they had "at least 
one penetrating encounter" (cf. Fig.7 in their paper noting 
that they are plotting the inverse of our distance measure). 
However, we no te that the definititio n for a satellite's virial 
radius used by ITormen et alJ lll99cf) differs to ours; they 
define the virial radius to be the satellite's radius at the 
moment it "merges" with the host halo, whereas we calculate 
the satellite's radius for each snapshot we have along its orbit 
within the host halo. This naturally leads to smaller radii as 
most of the satellites loses mass as it orbits within the denser 
environment of the host (c.f. definition for satellite radius in 
Section ^. 1,11 footnote 2), and as a result the distribution of 
relative distances peaks at larger separations. 

Relative velocities between satellites can also enhance 
the impact of interactions on mass loss - the slower the en- 
counter between a pair of satellites, the longer the timescale 
over which damage can be done. In Fig.|S]we show the distri- 
bution of relative velocities for pairs of satellites, normalised 
by the velocity dispersion of the host halo. As before, we 
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Figure 5. Distribution of the relative velocity of two satellites 
normalized by the host's velocity dispersion; for each host halo 
we stack the data from all available outputs. 



Figure 6. Measured fractional mass loss in between two consec- 
utive outputs versus the predicted mass loss as given by equa- 
tion 1121 . The solid line represents a X 1 / 3 power law whereas the 
dashed line indicates a simple 1:1 relation a; 1 ' . 



Table 2. Best fit parameters for relative distance distribution. 



Halo 


X 


O"0 


# 1 


0.097 


0.536 


# 2 


0.111 


0.510 


# 3 


0.158 


0.567 


# 4 


0.130 


0.550 


# 5 


0.104 


0.568 


# 6 


0.137 


0.564 


# 7 


0.147 


0.543 


# 8 


0.119 


0.540 



Table 3. Best fit parameters for relative velocity distribution. 



Halo 


w 


O"0 


# 1 


1.390 


0.563 


# 2 


1.319 


0.544 


# 3 


1.385 


0.572 


# 4 


1.323 


0.548 


# 5 


1.386 


0.564 


# 6 


1.358 


0.548 


# 7 


1.396 


0.562 


# 8 


1.292 


0.534 



stack data for all available outputs for each system, but we 
now fit the distributions with a Gaussian; 



n(w) — 



s/2r 



: exp 



Ulof 



(2a* 



(7) 



Here w = V rc \/cr^° st is the relative velocity of two satellites in 
terms of the velocity dispersion of the host halo. This figure 
suggests that there is no correlation of peak value with age, 
in good agreement with the best-fit parameters presented in 
Table El 

In summary, our analysis indicates that slow and/or 
close penetrating encounters between pairs of satellite galax- 
ies are relatively rare events. We have checked to ensure that 
our decision to stack all available outputs does not bias our 
result by masking a potentially interesting signal; however, 
we can confirm that the results are unaffected whether we 
construct the distribution from data obtained at a single red- 
shift (e.g. final redshift z — or formation redshift zf or m)- 



3.3 Mass Loss induced by Satellite-Satellite 
Interactions 

We have defined a physically motivated quantitative 
measure of interactions between satellite galaxies in the 
form of the IIM (equation EJ. Our detailed investigation of 
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IIM values in the previous sections provided great insight 
into the relevance of interactions in general; we were able 
to demonstrate that satellites encounter other satellites 
with the same relative velocities and separations at all 
times during the formation of a cluster, but the interaction 
measure IIM is higher in younger systems. However, it is 
difficult to conceive of a means to reconstruct IIM values 
for satellite galaxies from observational data. Whereas IIM 
values can readily be evaluated in TV-body simulations 
providing a gauge for the presence and importance of 
interactions, respectively, we prefer to construct a new 
measure for quantifying the impact of interactions more 
applicable to observational data sets. The mass loss suffered 
by a satellite as a result of the interactions would seem 
a promising approach; it can be probed observationally, 
such as in the field of "galactic archaeology" where tidally 
strip ped (stellar) stream s have proven to be a powerful tool 
(e.g. iHelmi et alJll999ft . However, understanding the evo- 
lution of satellite galaxies is complicated because changes 
are driven not only by the tidal field of the host (as shown 
by iKnebe et alJ 12005ft but also by more subtle processes 
such as the time evolution of the underlying host potential. 
Explicitly accounting for such time dependency gives better 
agreement with self-consistent modeling of satellites in 
the integrals-of-motion space, but there still remains a 
certain amount of disagreement betwe en the observ e d and 
measured mass losses; for example, IKnebe et al.1 ll2005ft 
speculated that this can be attributed to either the shape 
of the host and/or interactions with companion satellites. 
Using the ideas and prescriptions developed in Section [3. II 
we now extend our analysis to place constraints on the mass 
loss that can be induced by satellite-satellite interactions. 



A satellite i suffers a mass loss of 
AM' = M l (ti) - AT(t 2 ) 



(8) 



between two consecutive outputs ti and £2- We wish to relate 
a fraction of this mass loss to interactions between satellites 
and so we write AM 1 as the sum of the mass loss induced 
by interactions with other satellites, AM? at , and with the 
host halo, AM£ ost ; 



AM i = AACt + AML 



(9) 



In order to break the degeneracy between AA£? at and 
AMJ Mt , we assume that the satellites are point masses M 
with velocity v, and hence write their momentum change 
Ap as follows: 



Ap = vAM + MAv = FAt , 
which can be re-arranged to give 
AM = (F ■ vAt - MAv ■ v) /v 2 , 



(10) 



(11) 



where we (numerically) confirmed that on average 
(F • vAt) as 7 {MAv ■ v) and can therefore simplify the equa- 
tion for mass loss to read 



AM' x [ At 

v 2 (t m ) 



(12) 



where t m is the midpoint between two outputs calculated 
using equation J1J and a is a "tuning factor" accounting 
for the approximate nature of our approach. From Fig. 2] 
we conclude that a ~ 1/3 is the most appropriate value 



(represented by the solid line) and this is the value we adopt 
in the following analysis. 

We now use equation © and the scaling relation equa- 
tion \12l to compute the mass loss suffered as a result of 
satellite interactions and the influence of the host, respec- 
tively: 



AMI 



AMI 



AM 1 



1 + 



if(t m ) 



(tm)) 



(13) 



AM 1 



1 + 



where we have further assumed that equation 112H holds for 
both the force due to satellite-satellite interactions and the 
force induced by the host halo. 

Over the years, a number of sophisticated prescriptions 
for modeling tidally driven mass loss have been developed, 
largely within the context of th e evolution of globular clus- 
ters i n external tidal fields (e.g. ISpitzerlll958l : iGnedin et alJ 
Il999l and references therein) but also for understanding the 



disruption 


of satellite galaxies in 


cosmological dark matter 


halos (e.g. 


Pcnarrubia fc Benson! 


2005; Havashi et al.l 


2003; 



first order approximation for the mass loss, we will demon- 
strate that our formulae lead to qualitatively correct results 
and predictions with the right order of magnitude; a more 
thorough study and the development of a full theoretical 
model for mass loss in cosmological dark matter halos will 
be dealt with in a companion paper. In the present study 
we concentrate on quantifying the importance of interactions 
for mass loss and their importance for analytical modeling 
in galactic archaeology. 

In the following analysis we use the average fractional 
mass loss per Gyr for a given satellite i 



, dM xi 

( ) 

x Mdt' 



Snow 

-y 

TV? 



Ml - Ml 



Ml 



At 



(14) 



where JVJ is the number of outputs available for that par- 
ticular satellite between the time it enters the host and the 
present; the time interval At is calculated for two consecu- 
tive outputs. We are using equation I13H to split mass loss 
due to encounters and the influence of the host. The result- 
ing distributions for average mass loss per Gyr are shown in 
Fig. |7| This figure demonstrated that the mass loss induced 
by encounters between satellite galaxies can be as important 
as the tidal stripping of mass by the host potential in dynam- 
ically young systems. However, as the system becomes more 
relaxed, the relevance of such interactions becomes progres- 
sively less important and a significant fraction of the mass 
loss can be directly ascribed to the tides induced by the 
host. Table 2] accompanies Fig. [7] here we have calculated 
the mean of the average mass loss per Gyr for all satellites 
in a given host halo: 



dM 1 ^ dM i 

^Mdt' ~ l^^Mdt' 

From Table H we infer that the mass loss induced by 
satellite-satellite encounters can amount to as much as 45% 
of the total mass loss experienced by a single satellite. Even 



(15) 
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1.0 0.1 0.2 0.3 0.8 0.1 0.2 0.3 0.4 
average mass loss per Gyr average mass loss per Gyr 



Figure 7. Frequency distribution of the average (fractional) mass 
loss per Gyr. The thick solid line shows mass loss due to satellite 
interactions and the dashed line due to the influence of the host 
alone. 



Table 4. The mean average fractional mass loss per Gyr when av- 
eraging over all times and all satellites in a halo. The last column 
measures the contribution from satellite-satellite interactions, i.e. 



/ dM \ a dM \ 

\TZd~i ' sat ' \ Mdt ' total 



halo 


/ dM \ 
\Mdt' total 


/ dM \ 
\Mdt/ sat 


/ dM \ 
\Mdt' bost 


/sat 


# 1 


0.132 


0.030 


0.102 


0.23 


# 2 


0.123 


0.033 


0.090 


0.27 


#3 


0.143 


0.039 


0.104 


0.27 


#4 


0.126 


0.038 


0.088 


0.30 


# 5 


0.147 


0.050 


0.097 


0.34 


#6 


0.137 


0.041 


0.096 


0.30 


# 7 


0.141 


0.054 


0.087 


0.38 


# 8 


0.123 


0.055 


0.068 


0.45 



though the integral interaction measure IIM (as defined by 
equation 0) and its distribution in Fig.[2]indicated a rather 
low importance of such interactions the conversion to mass 
loss reveals a more pronounced influence due to the observed 
power law scaling AM oc (FAt) a with a ~ 1/3. However, 
the results are robust to changes in the power-law index, e.g. 
changing the exponent from 1/3 to unity gives us the range 
from 15% mass loss due to interactions for the oldest host 
(halo #1) up to 40% for the youngest system (halo #8). 



3.4 A Test Scenario: Host Halo #8 

Two questions remain unanswered: 

(i) Why do we observe a higher mass loss due to interac- 
tions in younger systems? 

(ii) Can our results be verified? 

The most natural approach to addressing these questions 
involves explicitly tracking the mass loss of an individual 
satellite as a function of time and factoring out the influ- 
ence of the other satellites. To do this, we have performed 
two additional simulation runs of halo #8, both starting at 
its formation redshift z = 0.3. In the first, we have removed 
all halos bar the progenitor of the z = host halo and one 
particular satellite that happened to have a rather high in- 
teraction value of IIM = 8.3; our analysis indicates that 
about 40% of its average mass loss per Gyr was induced by 
interactions with other satellites. In the second run, we re- 
moved only those subhaloes that have not merged with the 
host's progenitor at redshift z — except our "test satel- 
lite". We refer to these runs as "fully cleaned" (the former, 
including only the host and our test satellite) and "cleaned" 
(the latter, also including the massive subhaloes merging 
with the host) respectively. 

A visual impression of the initial setup of the cleaned 
run is given in (the upper left panel of) Fig. HUI which nicely 
demonstrate the "smoothness" of the cleaned simulation. 
The satellite in question is marked by a blue circle. 

For each of the three runs we closely follow the mass loss 
history of this satellite and the resulting curve (normalized 
to the initial mass) is presented in Fig.|S] Note that we iden- 
tify the set of particles that are bound to the satellite at the 
initial time and we explicitly track these particles through 
subsequent snapshots, checking what fraction are bound to 
the satellite at later times. This avoids any difficulties that 
may arise from attempting to identify the bound mass of the 
satellite as it is identified in subhalo catalogues constructed 
for consecutive snapshots. For each available output we find 
the new satellite centre by using the centre-of-mass of the 
innermost bound particles as a first guess for the central 
density peak. We then iteratively remove all of the satel- 
lite's particles that are not gravitationally bound. Fig. |H] 
shows that the mass loss suffered by the satellite is signif- 
icantly reduced in both of the cleaned runs. However, we 
stress that removing all substructure not only affected the 
satellite under investigation but also the overall dynamics of 
the host, especially for the fully cleaned run: host halo #8 
can be classified as a violent (triple) merger in the original 
cosmological simulation and removing all of its progenitors 
must clearly leave an imprint on its (internal) dynamics. 
Nevertheless, we note that although the "cleaned" run re- 
tained most of its high mass substructures, we observed a 
trend for the mass loss to decrease. 

This simple test has demonstrated that the principal 
driver for interaction induced mass loss is the mass spectrum 
of the substructure halos. We may strengthen this statement 
further by considering the cumulative circular velocity dis- 
tribution presented in Fig. [5] Noting that the maximum cir- 
cular velocity is a reasonable measure of halo mass, we con- 
clude from this figure that a sizeable fraction of the mass in 
satellites in each of the eight host haloes is bound to high 
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Figure 10. Several snapshots of the actual simulation of halo # 8 (right panel) and the re-simulation cleaned of all (but one) subhalo 
not ending up in the host at z = 0. The subhalo visible to the lower right of the host happens to be a foreground objects not interfering 
with the system under investigation. 



mass systems roughly half way through their evolution (i.e. 

t — (tform + to)/2). 

These massive subhaloes are likely to be responsible for 
the increased mass loss due to interactions, as already in- 
dicated by Fig. |H| We can investigate this assertion by de- 
termining the fraction of the IIM measure that is due to 
massive satellites. As already mentioned in Section 13.11 we 
can reproduce the distributions of IIM values presented in 
Fig.|5|by including only those satellite galaxies that are more 



massive than 1% of the host's virial mass 4 . This indicates 
quite clearly that the IIM values are dominated by more 
massive systems and the contribution of low-mass satellites 
is negligible. 

4 Using an empirically derived scaling relation between satellite 
mass and maximum circular velocity, i.e. v c [ IC oc A/ 1 / 3 , the mass 
limit of 1% M vir . 

host corresponds to a cut in circular velocity at 
around 20% u c irc,host (cf- Fig.[5J. 
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0.15 0.20 
redshift z 



Figure 8. Tracking the mass of the subhaloes indicated by the 
red circle in Fig, 1101 from initial redshift z = 0.31 to z = 0.08 
through the actual ("normal") simulation as well as out two test 
cases described in Section 13.41 



100 




Figure 9. Cumulative circular velocity distribution for all sub- 
haloes in the eight host halos half way through the evolution from 
redshift z = z toIIn to z = 0. The solid line represent the power 
law x~ 2 . 



4 CONCLUSIONS 

The hierarchical manner in which structure in our Universe 
forms - from the bottom up, through mergers and accre- 
tions - implies that interactions between galaxies (and con- 
sequently, between their dark matter haloes) are common- 
place. These interactions have been invoked to explain , for 
example, galaxy transformation (e.g. lMoore et al.ll99ctl. the 
excha nge of angular momentum (e.g^T^^^^^^^tothio^ 
1987), the triggering of star bursts (e.g. Mihos & Herncmist 
W94) and morphological change JsteuT^net^^NErTOrrol 
2002). The impact of interactions between satellite galax- 
ies and their massive host on the structur e of the satel- 
lites has been studied in some detail (e.g. lHavashi et alJ 



2003), but less well understood is the role played by interac- 
tions between satellite galaxies; in other words, the impact 
of satellite-satellite interactions. This important topic has 
formed the basis of this paper. 

We have denned an integral interaction measure (IIM) 
that allows us to quantitatively measure the importance of 
interactions between satellite galaxies for their mass loss. 
Our definition allows us to gauge the relative contributions 
of the host potential and other satellites for the mass loss 
suffered by an individual satellite. We have shown that the 
distribution of IIMs for a population of satellites within a 
cluster mass dark matter halo can be characterised as log- 
normal, and that the peak value (or mode) correlates with 
the age of the host system - typically the younger the host, 
the larger the peak IIM. Moreover, we note that the relative 
width of the distribution is broader in younger systems. We 
were able to confirm that the most significant contribution 
to the interaction measure comes from massive companion 
satellites which naturally explain the correlation with host 
age: subhaloes in young clusters have larger masses relative 
to the host since they have not been tidally disrupted yet 
which is validated by the observation that our younger hosts 
have a higher mass fraction in satellites. However, the IIM 
values are generally much less than unity, implying that the 
bulk of the mass loss suffered by a satellite is driven by its in- 
teraction with the host potential. We have also shown that, 
in those cases where the IIM is large, it cannot be due to 
single encounters; rather, it is built up through a series of 
many encounters. 

Our inv e stigat ions have also extended the result of 
iKnebe et al] (|2004Tl by demonstrating that not only are 
penetrating encounters between satellite galaxies relatively 
rare events over the "lifetime" of a cluster 5 , but that the 
timescale of such encounters is short, i.e. the relative veloc- 
ities are typically of order the ID velocity dispersion of the 
host. This result may be of interest to those engaged in de- 
veloping semi-analytic models of galaxy formation because 
we might expect the severity of encounters between satel- 
lites to be important for the efficiency of starbursts arising 
from tidal interactions. 

Finally we have proposed a simple empirical model for 
separating the respective contributions of the host poten- 
tial and interactions with other satellites for the mass loss 
suffered by a satellite. Our model suggests that mass loss 
driven by satellite interactions can be significant - in the 
particular test case we considered, we have shown that a 
given satellite can lose as much as ~ 40% of its initial mass. 
This may appear surprising at first, but we have shown that 
the IIM is a cumulative measure and so while damaging en- 
counters are relatively rare occurrences, a large number of 
"weak" interactions can affect the structure of a satellite and 
drive the mass loss it suffers. However, we stress that this 
empirical model should be taken as simple guide providing 
an "order-of-magnitude" estimate of the mass loss, and a 
more sophisticated model is required; this will be the focus 
of future work. 



5 i.e. since its formation redshift. 
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